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tical properties of the black hole degeneracy spectrum in loop quantum gravity. In 
particular we will study the persistence of the observed effective quantization of 
the entropy as a function of the horizon area. We will show that this quantization 
disappears as the area increases despite the existence of black hole configurations 
with a large degeneracy. The methods that we describe here can be adapted to the 
study of the statistical properties of the black hole degeneracy spectrum for all the 
existing proposals to define black hole entropy in loop quantum gravity. 
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I. INTRODUCTION 



The study of black hole (BH) entropy within the framework provided by loop quan- 
tum gravity (LQG) is an interesting issue that illuminates important aspects of quantum 
gravity. The modeling of black holes by using space-times admitting isolated horizons as 
inner boundaries, and the subsequent quantization of this sector of general relativity, has 
been extensively explained in the literature 111441. The resulting description provides a 
clear identification of the quantum BH degrees of freedom so that the standard quantum 
statistical definition of the entropy can be used. 

For small black holes the detailed behavior of the entropy as a function of the horizon 
area has been explored in BSHZl. A striking observation made in these papers is the fact 
that, in addition to the expected linear growth, the entropy displays a distinct staircase 
structure that amounts to its effective quantization. This is surprising because the spec- 
trum of the area operator is not equally spaced. A detailed study of this phenomenon has 
been undertaken by resorting to combinatorial methods -in particular the use of gener- 
ating functions- and number theoretic ideas. These have been described in ||7l- [T0ll . 

The so called black hole degeneracy spectrum is a way to encode the detailed infor- 
mation about BH configurations and their contributions to the entropy. In effect, the 
entropy can be computed as the integral of the black hole degeneracy distribution H5l[T0|. 
When this picture is used, the results on the entropy quantization manifest themselves 
as a distinct peak structure in the degeneracy spectrum (see Fig. [T]). This fact led to the 
identification in reference [8] of a peak counter -a function of the punctures of the spin 
network describing a BH state at the horizon- that efficiently labels the configurations 
contributing to a given peak. An alternative way to do this has been given in |[T0| as 
well as a generating function that singles out peak configurations. The main goal of the 
present paper is to use this master generating function to derive some important statistical 
information about the peaks in the degeneracy spectrum and discuss its physical impli- 
cations. The reason why we follow a statistical approach is the fact that an inspection of 
the nature of the degeneracy spectrum shows a combination of a simple coarse grained 
structure and a complicated detailed behavior as can be seen in Fig. [l] 

An important feature of our approach is the use of very strong results in Combinatorics 
that show a particular type of convergence to a Gaussian model when certain subsets of 
BH configurations are chosen. The obtention of the relevant statistical parameters, the 
mean and the variance, can be efficiently done in terms of the above mentioned gener- 
ating functions. We want to point out that some of the methods that we will use in the 
paper are particular applications of general theorems in Combinatorics (see the excellent 
book by Flajolet and Sedgewick [llj) suggesting that some features in the behavior of 
the entropy -in particular its effective quantization for small areas- are actually of a very 
generic nature. This is also supported by the fact that this phenomenon has been seen 
in all the different proposals found in the literature and is insensitive to the implementa- 
tion of the so called projection constraint [lOJ. We want to mention here that we restrict 
our analysis in the main body of the paper to the prescription given by Domagala and 
Lewandowski in ||T2|| to compute the entropy. In any case, we will show that our meth- 
ods can be easily adapted to deal with the other countings appearing in the literature, in 
particular the SU(2) proposal of [|3>|4J (see also lil3llT4l '), and our conclusions extended to 
these cases. 

The statistical information that we obtain has a direct physical application. By smooth- 




FIG. 1. Plot of the black hole degeneracy spectrum D (a), in units of 1 0^^, for a range of area values 
(in units of 47Ty£p, where y is the Immirzi parameter and l-p the Planck length). The band structure 
can be traced all the way back to the smaller values of the area. 



ing out the peaks in the degeneracy spectrum (or, rather, the steps in the entropy) and 
describing them as Gaussians (or better by Gaussian distributions, that can be written in 
terms of the error function erf) it is possible to obtain a smoothed representation for the 
black hole entropy The low area behavior, that has been studied so far in the literature, 
is captured in a very effective way by this model. It is possible to show that the interest- 
ing structure of the entropy seen for small black holes disappears in an area regime for 
which the smooth approximation is still valid. However, the preceding analysis does not 
exclude a revival of the entropy quantization for larger areas (or in the asymptotic limit) 
because the smoothed model fails to reproduce the exact value of the Immirzi parameter 
y and, hence, the correct growth rate of the entropy (although very good approximations 
for y are obtained in practice). 

The lay out of the paper is the following. After this introduction we will devote Sec- 
tion |ll] to give the basic definitions related to the entropy and the black hole degeneracy 
spectrum. In Section III we will study the statistical properties of the peaks by introduc- 
ing their moment-generating function. We will obtain the mean and the variance for the 
peak distribution and discuss the computation of higher moments. The approximation 
obtained by modeling the steps as Heaviside step functions (with discontinuities located 
at the area values given by the mean value of the areas associated with the peaks) will 
be discussed next in Section IV As we will see, the approximation obtained in this way 
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reproduces the behavior of the entropy for small areas remarkably well, though it is not 
suitable to understand the origin of the staircase structure itself. This can be better done 
by taking into account not only the mean but also the variance. We devote Section |V| to 
this issue. As we will see the low-area structure of the entropy can be neatly understood 
in this setting. Furthermore, we can also explain in quantitative terms how this structure 
disappears when the area increases. 

The basic approach discussed in the first part of the paper can be improved in several 
ways. One of them consists in further partitioning the space of black hole configurations 
by introducing extra peak counters. A particularly simple description can be found by 
using two of them. This will allow us to explain, at least in the limit of small areas, the 
appearance of discrete substructures in the peaks of the BH degeneracy spectrum. The 
details of this are described in Section VI We also give there a quantitative comparison 
of the peak counter found in [8] with other possible choices and conclude that it is the 
best one. We end the paper in Section VII with our conclusions and some details relevant 
for the extension of our methods to the SU(2) formulation of ||3l. A number of technical 
issues are left for the appendices. If not stated otherwise, areas in the paper will be given 
in units of Anyl^. 



II. BLACK HOLE ENTROPY: BASIC DEFINITIONS 

As we have mentioned in the introduction, some details in the entropy behavior as 
a function of the area are insensitive to the counting scheme that one chooses to follow 
(within the family of LQG inspired models). For the sake of concreteness most of the com- 
putations and results presented in the paper correspond to the Domagala-Lewandowski 
(DL) implementation [|T2| of the original proposal of Ashtekar, Baez, Corichi and Krasnov 
||ll|2|. However we will briefly discuss at the end of the paper the relevance of our results 
for the recent SU(2) proposal of [i3n4j. 

An extended discussion of the number-theoretic and combinatorial methods that we 
will employ here can be found in |fTO|, in particular the notation and definitions that we 
use in the paper. Nonetheless, and for the benefit of the reader, we give here the basic 
definitions that will be used. 

In the DL approach, the entropy S(a) (respectively S*(a), when the so called -projection 
constraint is ignored) of a quantum horizon of classical area a is given by 

S(a) =log(1 +aT(a)), (repectively S,(a) = log(1 + aT,(a))), 

where 9^(a) (respectively 91* (a)) is the number of all the finite, arbitrarily long, sequences 
(mi , . . . , trn ) of non-zero half integers, such that: 

N N N 

^ TUi = 0, lY_ VimiKlmil + l) < a (respectively 2 Y_ VlmiKlmil + 1) < a). 
1=1 1=1 1=1 

The condition rai = is known as the projection constraint. The computation of 
both S(a) and S*(a) can be efficiently performed in terms of the sets C(a'), a' < a, of 
the allowed configurations for each area a' — Y-i^i\/Vi ^ spffi'^'^'^) belonging to the 
spectrum of the LQG area operator. Here, as pointed out in [7J, qi G N U {0} and pt are 
the square-free integers (pi = 2, p2 = 3, pa = 5, p4 = 6, etc.). A configuration c e C[a), 
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as defined in IITOll , is a (finite) multiset c = {(k, Ni<)}in which each integer k G N appears 
N]^ times (with Ni^ 6 N U {0}). The set C of all possible BH configurations is defined as the 
union of the configurations corresponding to the different area values 

e:= U e(a). 

There are several functions f : C ^ M defined on the space of configurations with a 
clear physical interpretation that we will use extensively in the following: 

• For any configuration c = {(k, Nj^)}, N (c) = represents the number of punc- 
tures defined by a spin network piercing the horizon, K(c) = ^j^kN]^ is (twice) 
the total spin, and A(c) = 7k[k + 2)Nk the area of the BH induced by c. We 
will use also some other functions defined on C such as the "peak counter" P(c) := 
3K(c) + 2N (c) . These functions satisfy the bound P/3 < A < 2P/5 or, equivalently, 
5A/2 < P < 3A. 

• As explained in pTOl, the degeneracy d(c) of a configuration c G C(a) allows 
us to compute the number D(a) = Y.cee{a) ^^^"^ arbitrarily long sequences 
(mi , . . . , raN ) of non-zero half integers, such that: 

N N 

}^mi = 0, 2}^ V|miK|mi| + 1) = a. 
1=1 1=1 

In terms of the so called BH degeneracy spectrum D ( a), the BH entropy is given by 

expS(a) = 1 + D(a') . 

a'<Q 

• When the projection constraint is ignored, the degeneracies d* (c) of a configuration 
c G C(a) give us the number D*(a) = Y.cee(a) arbitrarily long sequences 
(mi , . . . , mN ) of non-zero half integers, such that: 

N 

2}^ V|mi|(|mi| + 1) = a. 
1=1 

The BH degeneracy spectrum D*(a) in the case when the projection constraint is 
not considered can be used to compute the value of the entropy 

expS,(a) = 1 + D,(a') . 

a'<Q 

For a given area a = + c\2y/V2 + ■ ■ ■ , the values of D(a) and D*(a) can be encoded 

in the coefficients of the generating functions G (z; xi , X2, . . . ) and G* (xi , X2, . . . ): 

D(a) = [z°][x?'x5^---]G(z;xi,X2,...) 

(oo oo 
1 -^^(z^- + z-'^Mx^^ 
i=1 n=l 

D*(a) = [x^'x^^ • • •]G*(xi,X2, . . . ) 

-1 



= [x?'xf 



i=l n=l 
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For each square-free pt, the terms (z^^ + z^^^jx^^^ and x^"^ appearing in the corresponding 
generating functions are built from the solutions {(k{:^, Unjln to the Pell equations (k+ 1 — 
Pi^y^ = 1 (see \\lOi for details). Here [z°][x^^x5^ ■ ■ ■]G(z;xi,X2, . . . ) denotes the coefficient 
of the z°x^' x^^ ■ ■ ■ term in a Laurent expansion of G (z; Xi , xi, . . . ) about z = 0, Xi = 0, . . . 

It is important to notice that the family {C(a) C C : a G sp^Q^'^'^)} provides us with a 
partition of the configuration space C defined in terms of the level sets of the area func- 
tion C(a) = (a). If we are given any other function P in the configuration space (in 
particular P(c) = 3K(c) + 2N(c)) it is possible to define a different partition C = Up 
using the level sets CPp = P"^ (p). This means that the sets Cp(a) = P^^ (p) n A^^ (a) define 
a finer partition than either |J^ C(a) or IJ^ 7^: 

V a 

Notice that C(a) = [j^ Cp(a) and CPp = |J^ Cp(a). This fact can be used to compute the 
entropy as 

exp S(a) — 1 + ^ ^ D(a' I p) , exp S,(a) = 1 + ^ ^ D,(a' | p) 

p a'<a p a'<a 

where 

D(a|p)= Y_ d^c), D,(a|p)= ^ d,(c). 

ceCpta) ceep(a) 

This is so because 

expS(a) = 1 + D(a') = 1 + Y_ d(c) = 1 + }^ Y_ ^^^^ 

a'<a a'<acee[Q') Q'<a p ceCpla') 

= ^+^Y. L d(c) = l+}^}^D(a'|p) (II.l) 

p a'<Qceep(a'] p Q'<a 

and equivalently for exp S^, (a). Finally, it is important to notice that, when the partition is 
defined by the functions of the type P(cx, (3) := ocK + (3N (a generalized "linear" counter 
with positive integer coefficients), the numbers D(a | p) and D^,(a | p) associated with an 
area a = (]^^/p^ + ^2y/V2 H can be derived as 



D(a|p) = [z°][x?'xf ■■■][a^^]G(v,z;xi,X2,...), D,(a|p) = [x^^x^^^ ■ ■ ■][^p]G(v;xi,X2, 
from the master BH generating functions 



oo oo 



-1 



G[y,z;xuX2,...) := i^ -^^[z^^ +z-^^)y°'^^^^xf\ , (II.2) 

\ i=l n=l / 

/ oo oo \ 

G.(a/;xi,X2,...):= 1-2}^}^v«^n+Pxf^ . (II.3) 



1=1 n=l 



Notice that these generating functions are normalized in such a way that D (0) = D (0 1 0) 
DJO) =D,(0|0) = l,andD[0|p) = D,(0|p) =Oforp^O. 
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III. STATISTICAL PROPERTIES OF THE PEAKS 

The starting point of our analysis is to introduce a convenient partition of the space of 
black hole configurations that is adapted to the description of the peak structure seen in 
Fig. [Tjfor the BH degeneracy spectrum (or, alternatively, to the steps of the entropy). This 
partition is performed by introducing the peak counter P = 3K+2N defined above. A peak 
in the space of BH configurations CPp C C is defined as consisting of those configurations 
c corresponding to a pre-selected value p of P(c). We have then 

V P n 

For a fixed value of p there are, of course, configurations corresponding to different val- 
ues of the area (within a bounded range P/3 < A < 2P/5) and different degeneracies d(c) 
(or d*(c)); this is so because Tp = Cp(a). In fact, if one plots D(a,p) = ^ceep(a) '^f'^) 
or D*(a, p) = Y.cee^{a] d*(c) as a function of a (for a fixed value of p) one gets a regular 
structure with the form of peak in the black hole degeneracy, as shown in Fig. |2] It is 
obviously possible to reconstruct the degeneracy spectra that have already appeared in 
the literature |fl0| by adding up the contributions of these peaks for all the values of p. 



The regular shape seen in Figure |2] strongly suggests that a Gaussian approximation can 
provide a good description of the peaks. This fact leads naturally to the consideration of 
statistical methods to study BH entropy. In fact, this is the main theme of this paper. We 
want to emphasize from the start that we do not merely compute statistical parameters 
(the mean, the variance and, eventually, higher moments) by fitting Gaussian profiles to 
the peak data but, rather, obtain them exactly from the BH generating functions. In other 
words, we will not just use descriptive statistics, but employ the very powerful analytical 
tools available for a wide class of combinatorial problems (involving generating functions 
of the same type as the ones that we use in this paper [llj). This will allow us to make 
predictions regarding the statistical parameters of arbitrary peaks and use them to study 
the behavior of the BH entropy. 

A statistical treatment requires us to give a weight to each configuration. In our prob- 
lem this is naturally provided by the degeneracy d(c) -or, respectively, d*(c). The rele- 
vant objects to be computed are the expectation values of the powers of the area (taken 
as a random variable) conditioned by a fixed value of p: 

^^gy d(c)A'^(c) ^^gy d,(c)A^(c) 

E[A-|p]:= ^^^"^ , EjA-|p]:=^^f^— — — . (III.l) 

In the first case (where the projection constraint is taken into account) only even values 
of p have be considered because Y.ceTp ^^^"^ otherwise zero. As we will explain later, 
we will use the relevant moments defined by this formula to build a smooth approxima- 
tion for the shape of each step in the entropy. This will require us to "de-normalize" the 
distribution by multiplying it by the total peak degeneracy Y.c&p ^^^"^ X^ceTp 
The standard way to compute E [A"^ | p] and E* [A"^ | p] relies on the use of the so called 
moment-generating function associated with the random variable A. A remarkable feature 
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FIG. 2. Peak corresponding to p = 150 in units of lO^'^. By plotting the values of D*(a|p) one gets, 
for the largest degeneracy values, a roughly symmetric shape. The values corresponding to the 
largest degeneracies can be readily seen in the plot. Notice, however, that the peak extends to the 
left much farther than it does to the right. This phenomenon can be easily seen in a plot of the 
logarithm of D* ( a | p] that also displays a distinct subpeak structure to be discussed later (see Fig. 
Tol. 



of the combinatorial approach that we follow to study black hole entropy in LQG is the 
fact that this moment-generating function can be easily derived from the master gener- 
ating functions ( II.2[ > or (11.3) given above. We will start by looking at the case where the 



projection constraint is ignored. The incorporation of the projection constraint will be 
discussed afterwards. Though this problem is more complicated, there are no important 
conceptual differences as far as our treatment is concerned. 



III.l. Moment-generating function: Ignoring the projection constraint 

Let us take as the starting point the master generating function G* ( v; xi , . . . , ), defined 
in (II.3I, where the variable xt refers to the square free integer pt. By substituting Xi = 
g-VPi/as is standard in this setting IITOll , we obtain 

G,(v, s) := G,[y; e'^^, e'^^^ . . .) = ^ ^= • (III.2) 

1 - 1 v3^+2e-V'^f'^+2' 
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By construction, it is obvious IITOll that 



-sA{c) 



and, hence. 



EJexp(-sA) Ip] 



(in.3) 



Modulo normalizing factors, and the exchange s i-^ —s, the function g(s |p) := [a'P]G,,(v, s) 
is the standard moment-generating function used in Mathematical Statistics and, hence, 
logg(s Ip) is the cumulant-generating function used in Statistical Physics. Notice that 
our sign convention originates in the use of Laplace transforms to write down closed ex- 
pressions for the black hole entropy | IT5llT6| . By computing the derivatives of (III. 2) with 
respect to s at s = we can easily find all the expectation values for arbitrary powers of 
the area: 

EjA"ip] = (-ir — 



s=0 



In particular, the mean and the variance 



|x^p = EjA|p], = EjA^lp] -E^[A|p] 

of the area distribution conditioned by P = p, can be obtained in a straightforward way. 
Exact expressions (as closed functions of p) for [x^p , and the normalization factor 



can be found in Appendix |Aj In the asymptotic regime p ^ oo these objects follow very 
simple laws: 



~ [J, . p = (0.34959022 ■ 
ol^^a^.^ = [o.oooomi- 



)-P, 



(111.4) 

(III.5) 
(III.6) 



where = (0.77039825 ■ ■ ■ ) will denote the single real root the polynomial 2y^ + y^ — 
A closed expression for yo in terms of hypergeometric functions is given in [|T0| (we will 
discuss some details concerning this issue in Section VI). As we show in Appendix A 
the coefficients yi and appearing in (III.5I and (III.6) can be written in terms of yQ. The 
linear growth of the variance and the fact that the spacing between successive steps tends 
to a constant value strongly suggests that the steps will fade as the area increases. This 
will be shown in detail in Section fVl 

We want to mention here that the leading behavior of |j.*p ~ (0.34959022 ■ • ■ ) ■ P exactly 
coincides with the one obtained in [lOJ by using a completely different approach relying 
on a continuum approximation. We will provide an alternative proof of these asymptotic 
results in Section Ivl 
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III.2. Moment-generating function: The Domagala-Lewandowski approach 

When the projection constraint is considered, the starting point is the master generat- 
ing function G ( v; z, Xi , . . . , ) defined in (11.2 1. By substituting Xi — e^^v^ we obtain 

G ( V, s; z) := G (v, z; e^^^^, . . .) = ^ , . (in.7) 

The function G (a^, s; z) satisfies 

[z°][a'P]G(a/,s;z) = d(c)e-^^'^' 

ce^p 

and, hence, we have the following expression for the expectation value 

^^^"P^-^-^^'^^= [zO][.p]G(.,0;z) - ^"^-'^ 

By computing the derivatives of ( |IIL7 | with respect to s, at s = 0, we can easily find all 
the expectation values for arbitrary powers of the area 



E[A-ip] = (-ir- 





( 9^ 






G(v,s;z) ) 




s=0 / 



[zO][yp]G(^,0;z) 
In this case, the mean and the variance 

^p = E[Aip], = E[A2|p]-E2[A|p] 

of the area distribution conditioned by P = p can be obtained with some extra work due 
to the presence of the z-variable in G(v, s; z). Expressions for |a.p ,crp , and the normaliza- 
tion factor 

(Xp:=[z°][A/T']G(v,0;z) = _^d(c) 

ceyp 

can be found in Appendix |Bj In particular, it is possible to prove that ap = for all odd 
values of p and hence only the even values of P(c) = p have to be considered. In this 
case, in the asymptotic regime p = 2q ^ oo we have 



^2q ~ 2|x ■ q = 2 ■ (0.34959022 ■ ■ ■ ) ■ q , (ULIO) 
(T^q ~ 2o^ ■ q = 2 ■ [id.mQ')^^ 7 ■ ■ ■ ) • q . (ULll) 

The statistical treatment given in this section suggests two approximate models for 
the behavior of the black hole entropy as a function of the horizon area. In the first one 
the steps in the entropy are approximated by Heaviside step functions with jumps of 
magnitude ocp (or, respectively a*p when the projection constraint is ignored) located at 
areas given by the mean values |j.p. The second, improved, model will use smoothed 
steps given by the (integrated) Gaussian distributions of mean |j.p and variance cTp with 
height Kp. We discuss them in the following sections. 
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IV. USING THE MEAN: THE STAIRCASE APPROXIMATION FOR THE ENTROPY 

A coarse approximation for the exponentiated entropy exp S*(a) and exp S(a) can be 
obtained by assigning the sum of all the degeneracies corresponding to each peak to a 
single step located at the mean area value. This can be done by employing Heaviside 
step functions (denoted by 9 in the following) in several slightly different ways (obtained 
by using the asymptotic approximations for |j.^,p, [i^, a*p and ocp): 

oo oo 

^ cx*p9(a- |j.=,p) , y~ ap9(a- \Lp) , 

p=0 p=0 
oo oo 

^ a*p9(a- |a.p) , y~ ap9(a- \xp) , 

p=0 p=0 



where the first column corresponds to the case without the projection constraint. The 
validity of these approximations for small areas is clearly seen in Fig. |3] 



By proceeding in this way the entropy will obviously display a staircase structure be- 
cause we are approximating it as a sum of sharp, (asymptotically) equally spaced, steps. 
This means that this simplified approach will not be suitable to address the persistence (or 
lack thereof) of the structure seen in the entropy for small areas in the asymptotic regime. 
However, it can be used to estimate the value of the Immirzi parameter as y Tstairs 
because it provides a simple expression for the growth of the entropy as a function of the 
area. In fact, in the case where the projection constraint is neglected, we easily find 

TTTstairs := = (0.74615268 ■ ■ ■ )• (IV.l) 



This value must be compared with the one obtained by Meissner in IITSl 



7ty = (0.74623179- ■■)■ (IV.2) 

Notice that ystairs has nothing to do with the value of the Immirzi parameter Ttyflux = 
log 3 derived in the context of LQG from the evenly spaced flux area operator used in 
[[T7| which can be interpreted in terms of the Schwarzschild quasinormal modes [flTl [18| 



Though one could argue that the value ystairs derived for the Immirzi parameter in this 
approximation is quite good, the fact that it predicts a growing behavior, different from 
the true one, means that the entropy and its staircase approximation will diverge linearly. 



A convenient way to derive (IV.l) is by using Laplace transform techniques. Let us 
discuss, in the first place, the staircase approximation without the projection constraint. 
To this end we consider 



2v^ °° 1 
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5 10 15 Q 



FIG. 3. The upper part of the figure shows the exact value of the entropy S* ( a) for areas below 1 8 
(in units of 47ty£p) in the case in which the projection constraint is ignored. The lower part repre- 
sents the staircase approximation to the entropy given by Sf ^'''^(a) = log f '^p0(o- ~ 1^ " P) ) • 
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whose Laplace transform L[exp Sf^"^, s] can be computed in closed form (as a function 
of the complex variable s) 



p=5 



L(exr> S"*"*^"' si = - H ^ V - = - H 



The pole Sgtairs = ~ log(^o) /m^ is responsible for the exponential growth of exp Sf^^''®(a), in 
the regime a ^ oo, given by ( IV.1[ >. Notice that, in addition to this real pole (and s = 0), 



there are infinitely many others of the form — log(A'o) 2k7ti/(x, k G Z \{0} that account 
for the steps in this approximation for the entropy |[T6|. 

When the projection constraint is taken into account the configurations c with odd 
values of P(c) = p have zero degeneracy and hence only even values of p = 2q have to 
be considered. The staircase approximation is then 



oo 



The Laplace transform L (exp 5^*^^''% s) is given by 



£(expS^'^^-,s) = - + r^W 7in 2ZlV iv.(e-^'^-+^°g-°') 

where Liy^ denotes the polylogarithm of order Vi. The singularities of Lii/j(e^^f^*+'°s^''') are 
branch cuts starting at the same straight line Re(s) = Sgtairs in the complex s-plane as the 
singularities found for the case without the projection constraint: — log(vo) /m- + km/ [i, 
k G Z. Notice that the spacing between these points is half the one obtained when the 
projection constraint is ignored. This means that the width of the steps doubles in this 
case. The effect of the branch cuts is to modify the asymptotic behavior of the entropy 
by the addition of the expected logarithmic corrections, however, the linear growth is the 



same as before and the inferred value of the Immirzi parameter is still given by (IV.l I. 

The failure to reproduce the exact value for the Immirzi parameter in this approxima- 
tion stems from the fact that, for a given value of the area a, the model neglects to take 
into account contributions coming from peaks with p beyond the largest one satisfying 
\xp < a. It also misses some contributions coming from lower values of p (at least in the 
asymptotic regime of large areas). 

V. USING THE MEAN AND THE VARIANCE: SMOOTHED GAUSSIAN 
APPROXIMATION FOR THE ENTROPY 

An improved model for the black hole entropy can be obtained by approximating the 
steps by Gaussian distributions with mean and variance given by ( III.5[ > and (III.6 1. This 



will take into account the fact that the steps become wider with increasing values of p (an 
effect that can be readily seen by plotting the exact values of the entropy for small black 
holes as functions of the area). This is obviously relevant to study whether the staircase 
structure is present in the asymptotic limit. At this point it is just appropriate to quote 
from page 611 of the book by Flajolet and Sedgewick |[TT]| 
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"Many applications, in various sciences as well as in combinatorics itself, require quantifying 
the behaviour of parameters of combinatorial structures. The corresponding problems are now 
of a multivariate nature, as one typically wants a way to estimate the number of objects in a 
combinatorial class having a fixed size and a given parameter value. Average-case analyses usually 
do not suffice, since it is often important to predict what is likely to be observed in simulations or 
on actual data that obey a given randomness model, in terms of possible deviations from the mean 
-this signifies that information on probability distributions is wanted. [...] Indeed, it is frequently 
observed that the histograms of the distribution of a combinatorial parameter (for varying size 
values) exhibit a common characteristic "shape", as the size of the random combinatorial structure 
tends to infinity. In this case, we say that there exists a limit law. " 

In our case we have a multivariate combinatorial problem where both the area and the 
peak parameter P play a significant role. Furthermore, we have that the distribution of 
one of the parameters (the area of the peaks) displays a characteristic shape as the peak 
counter grows towards infinity As we will show in this section the methods appearing 
in [fTTjl will allow us to gather important information about the behavior of the entropy 
as a function of the area. In particular, we will see that a Gaussian law -reminiscent of 
the Central Limit Theorem of probability theory- plays an important role in the analysis 
presented here. 



V.l. Gaussian law for the peaks 

The key idea -in the case where the projection constraint is neglectec^ is to use the- 
orem IX.9 (page 656) of |[TT| for the generating function G^f^^, s) given in ( |III.2 1. The 
theorem tells us that the mean and the variance for P = p can be easily obtained in terms 
of the "analytic" mear|^m(f*) and variance t)(f*) of a function f* as 



= m(f Jp + 0(1 ) = • p + 0(1 ) , (V.l) 



The function 



< = »(f.)p + On)^(f<«|-(^<|)^)p.O(,). (V.2) 

f*(s) := — — , f,(0) = l, 
is given in terms of y*[s) defined by y*[0) = and Q*(v*(s), s) = 0, where 

^ oo 

Q*(^,s) := r = ^-2yy'''+'exp[-s^m + 2)) . 

The implicit function theorem allows us to obtain a power series expansion in terms of 
the variable s with coefficients given by derivatives of Q*(a', s) evaluated at v = a^q arid 



^ The case when the projection constraint is taken into account can be handled by adapting theorem IX. 12 
of 111]. 

^ Notice that the minus sign in our definition of m(f * ) originates in our sign convention for the variable s 
appearing in our moment-generating functions. 
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s = 0. Explicitly 



=^0 



2q; 



1,0 



with 



q*i,j := 



The results given by (V.ll and (V.2| are exactly the same that we have found above in 
equations ( IIL5[ > and (IIL6). In any case this is a very efficient method to compute the 
numerical values of the mean and the variance of the peak distributions. The theorem, 
however, provides us with another very important convergence result: The random vari- 
able 



with (normalized) distribution function 



0" 



T^*p(^) = Prob*(Xp < x] 



*p 



ceypnXpi((-oo,x]) 



converges, pointwise, to a Gaussian distribution 



.x'<x Z-ceypHXp '(x'J 



lim f^Jx] 

p^oo / 



e ^ dt , 



with a 0(1 /y^) speed of convergence. In terms of the area this fact implies that we can 
write 



Q'<Qceep(a') 



a'<Q 



y~ d*(c) = a*pF=^p ( — ^) 



ceTpnAp'(o,Q] 

In practice this tells us that each smoothed step, given by the function 



a 1-^ 



is a good approximation (see Fig. |4]) to the actual shape of the graph of the function 

a ^ D^a' Ip) 



a'<Q 



appearing in the definition of the entropy (II.l ). This approximation improves as p grows 
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FIG. 4. Comparison between the exact shape of the step X.a'<a D*(a' 1 120) and the correspond- 
ing smoothed Gaussian approximation ( 1 + erf 



around the value |j.*i2o and has a width of 2 cx^, 120- 



Q-H-*120 



The vertical band is centered 



If one compares, instead, the graphs of the functions 



<x 



2710". 



exp 



a 



and a 1-^ D,,(a|p) , 



corresponding to the peaks in the degeneracy spectrum, the Gaussian shape does not 
correspond in any way to an "envelope" of the actual peak defined by the degeneracies 
(see Fig. [l]), although the maxima appear roughly for the same value of the area and 
the widths match reasonably. We want to point out that the parameters of the Gaussian 
approximation have been obtained a priori from the moment-generating function. So we 
are not fitting the "peak data" to a Gaussian but, rather, deriving the statistical properties 
of the distribution that they define by relying on an exact statistical analysis. 



V.2. Gaussian approximation for the entropy 



The idea now is to model the entropy as a sum of smoothed steps like the ones shown 
in Figure |4j However, it is very important to be aware of the fact that the convergence 
of the individual steps to their Gaussian approximations does not guarantee the conver- 
gence of their sum to the actual value of the entropy. Let us consider then approximations 
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for the exponentiated entropy exp S^, ( a) and exp S ( a) obtained by adding Gaussian steps. 
In the case where the projection constraint is ignored these are 




Similar expressions (to be discussed later) hold for the case in which the projection con- 
straint is incorporated. Figures|5]and[6]show a comparison between these approximations 
and the actual value of the entropy for the smallest areas. As it can be seen, the agree- 
ment is excellent. Not only the height of the steps is reproduced with high fidelity but 
also their progressive smoothing. In the case when the projection constraint is taken into 
account the staircase structure is more evident. 



In spite of this remarkable agreement we know, as we have learned in the preceding 
section, that the asymptotic growth of the Gaussian approximation may not be the exact 
one (i.e. the value of ygauss derived here may differ from the true one). In fact, this 
will be shown to be the case. In order to study the asymptotic growth of the Gaussian 
approximation to the entropy it suffices to consider 

and study the singularity structure of its Laplace transform written in terms of the com- 
plex variable s. This is given by 

.(expsr",s) ^ 1 + y^o^l ^ (1 -erf {^^^)) (V.3) 

Notice that despite the fact that the Laplace transforms of the individual steps are entire 
functions in the complex variable s, the analytic extension of their sum (restricted to the 
values of s for which it actually converges) may, and actually will, have singularities. 
By looking at the inversion formula for Laplace transforms it is easy to see that these 
singularities determine the asymptotic behavior of the original sum. This is discussed in 
detail in Appendix [Cj The first two terms in (V.3 1 have a very simple analytic structure 



because they just have a pole at s = 0. The second sum in (V.3) is more complicated and 
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12 - 



10 




FIG. 5. The upper part of the figure shows the exact value of the entropy S^,[a) for areas below 
1 8 (in units of 47ty£p) in the case in which the projection constraint is ignored. The lower part 
represents the Gaussian approximation Sf^"^^(a) = log ^1 + ^ ^1 + erf 
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FIG. 6. The upper part of the figure shows the exact value of the entropy S(a) for areas below 1 8 
(in units of 47ty£p) in the case in which the projection constraint is taken into account. The lower 

part represents the Gaussian approximation S8'^"®^(a) = log ^1 + ^ (^1 +erf ^^^^^^^ . 
The staircase structure decays more slowly now as discussed in the text. 
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may converge or diverge depending on the values of s. It is easy to see that it converges 
for all s G C such that |arg(s)| < 37r/4 (the range of the argument is taken to be arg(s) G 
(—71, n]). Inside the wedge |arg(s)| > 3n/4 the series converges for values of s to the right 
of the hyperbola (see Figure 11 in Appendix |C]) 



Reims' 



M-s - log Vo ) = . 



This divergence (in the region |arg[s)| > Sn/A) is due to a term of the form 



:iov^ + 3)s^v: 



p=5 



(V.4) 



(V.5) 



This means that by subtracting this expression from the series that we are looking at, 
we get another series that converges in the full wedge |arg(s)| > 37r/4 to a function h(s) 
(with no singularities). The sum (V.5 1 can actually be performed in closed form to get a 
meromorphic extension to |arg(s)| > 37r/4 of the function that it defines inside its region 
of convergence. This is given by 



2exp 



5|j.s 



^1 -exp (^ys^- M-s-logvo)^ s 



(V.6) 



v^(10v^ + 3) ( 1 



The analytic extension of the Laplace transform (V.3| to the wedge is then given by the 
sum of the first two terms in ( V.3 1, the function hfsjand (V.6 1. Hence the singularities of 
the Laplace transform ( V.3[> ar e s = and those of (V.6 1. These are isolated simple poles 
located on the hyperbola (V.4) defined above (see Fig. 11 1 and given by the condition 



— |a.s — log ^ro = Ikni , k G Z , 



The single real pole at 



s --^ n 

^gauss — ^2 ' 



1 + ^logvo 



dictates the asymptotic growth of the entropy in this Gaussian approximation and gives 
an improved estimate of the Immirzi parameter 



Tty gauss = Sgauss = (0.74623087 • 



(V7) 



to be compared with the actual value ny = (0.746231 79 ■ • • ). As it can be seen, the value 
of Ygauss is better than ystairs obtained in the preceding section by using the staircase ap- 
proximation but still not the true one y, as they differ starting at the sixth decimal figure. 

The decay of the staircase structure is dictated by the two poles with the smallest non- 
zero imaginary parts. These are given by 



■2 



0" 



1 + ^(logvo ±2m) 



[0.70084660 • • • ) T (1 7.97653845 
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The magnitude of the imaginary part in the previous expression is very close to In/ \i. = 
(1 7.97300040 ■ ■ ■ ) expected. Finally the comparison between Sgauss and the real part of 
tells us the decay rate of the staircase structure of the entropy. This given, essentially, 

by 

exp ( - (sgauss - Re(s,±)) a) = exp ( - (0.04538426 ■ ■ ■ )a) (V.8) 

which means that the steps should have faded significantly for areas with an order of 
magnitude given by 1 / (0.04538426 ■ ■ ■ ) ~ 25. We show in Fig. [7| the exact behavior of 
the entropy for two different area intervals in the case when the projection constraint 
is ignored and the corresponding Gaussian approximations. It can be readily seen both 
the accuracy of the Gaussian approximation in this regime and the decay of the staircase 
structure (essentially absent for areas around 50). Another interesting feature that can be 
seen in Fig. [7|is the fact that as the density of the area spectrum increases, the jumps in 
the values of the entropy for consecutive area eigenvalues become smaller and smaller. 
Hence the entropy is better and better described by a smooth curve (in fact a straight line). 
It should be pointed out that the BH degeneracy spectrum at this regime still shows a dis- 
tinct peak structure produced by configurations of very large degeneracy that, however, 
give an almost negligible contribution to the total degeneracy of the individual peaks for 
large areas. 



When the projection constraint is taken into account the Gaussian approximation is 



and the singularities of its Laplace transform -that again control the asymptotic behavior 
of the entropy- are encoded in the series 



that plays the same role as ( V.5[ > when the projection constraint was ignored. The singu- 
larities satisfy now 

— s — |j.s — log vo = km , kGZ. 

The single real branch point is located at Sgauss/ which is precisely the place where the real 
pole in the case without the projection constraint is placed. The change of the singularity 
nature (branch points instead of poles) means that the asymptotic behavior of the entropy 
in the Gaussian approximation (with the projection constraint) will acquire the expected 
logarithmic corrections. The decay of the staircase structure is controlled by the two 
singularities with the smallest non-zero imaginary parts. These are given by 



s± := ^ M - Wl + ^(logvoiTti) ) = (0.73488453 ■■■) T (8.98835516 ■■■ )i • 
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50.0 50.5 51.0 51.5 Q 52.0 



FIG. 7. Comparison between the exact values of the entropy S* and their Gaussian approximation 
for two different area intervals. The necessary computations are made possible by the introduc- 
tion of the master generating functions ( |II.3| . We can see that the Gaussian approximation works 
remarkably well in both cases. The fading of the oscillations in the entropy for the largest area 
interval is evident. Similar plots for area intervals around 70 show an essentially linear growth. 
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Two important points must be mentioned now. First we notice that the real part of these 
singularities is much closer to Sgauss than in the case without the projection constraint. 
The decay rate of the staircase structure is controlled now by 

exp ( - (sgauss - Re(s±)) a) = exp ( - (0.01134633 ■ ■ ■ )a) . (V.9) 

This means that the staircase structure of the entropy will be perceptible for areas roughly 
four times as large as those corresponding to the non-projection constraint case. Hence, 
similar results to those shown in Fig. |7] are obtained with the projection constraint for 
areas around 200 in our units. Finally, the imaginary part of the singularities -that con- 
trols the width of the steps- is essentially half the value that it had when the projection 
constraint was ignored and hence the steps are twice as wide. 

VI. OTHER PARTITIONS 

An interesting question that naturally arises concerns the comparison of the results 
obtained with the peak counter P = 3K + 2N that we have used throughout the paper 
with the ones obtained with other partitions defined by linear counters of the type 

P(a, (3) := aK+ |3N , a, |3 G N. 

An important point that we want to emphasize here is the fact that all these counters 
provide partitions of the space of black hole configurations, and hence the full black hole 
spectrum can be recovered by using any of them (and taking into account all its possible 
values). However, some of them may be better suited to understand or isolate specific 
features of the spectrum (such us the observed bands for low areas). 

By reasoning as in [8], and relying on the solutions to the Pell equation, it is actually 
possible to find other counters, such as P(5, 4) = 5K + 4N or P(7, 5) = 7K + 5N, that could 
be potentially useful to understand the degeneracy spectrum. In particular, we have seen 
that P(3,2) works very well for small areas but leads to a Gaussian approximation that 
underestimates the value of the Immirzi parameter. It is then natural to wonder if a 
better counter could exist that provides a better estimate for y and still explains the low 
area behavior of the entropy. Even if the low area behavior is not captured by such a 
counter it could be used to understand the asymptotic limit of large areas. In this section 
we compare the P(3, 2) counter with P(cx, (3) and show that our original choice is the best 
one. This does not mean that these other counters are not useful themselves. In fact, we 
will show that they can also be used to refine the partition provided by P(3, 2) and study 
the interesting substructure of the peaks defined with the help of P(3, 2). 

VI.l. Assessing the "goodness" of P(3, 2) = 3K + 2N 

The level sets Tp(cc, |3) = P(cx, |3)-^ (p) of the function P(cc, |3) = aK + |3N lead to the 
partition C = Up ^p('^> of the configuration space. It is important to point out that 
the partitions defined by P(na, n|3) are equivalent to the partition defined by P(a, |3) for 
any positive integer n. This means that we can consider only values of a and |3 that are 
coprime, i.e. such that gcd(a, |3) = 1 . 
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By proceeding as in previous sections we arrive at the following peak generating func- 
tior£| 

G, [y, s; a, (3) := ^ == . (VI.l) 

The mean and variance of the area distribution conditioned by P((x, |3) = p can be easily 
derived from this generating function and have an asymptotic behavior given by 

lXp((x, |3) ~ M,(a, |3) ■ p , o-J(a, |3) ~ 0^(a, (3) ■ p , 

where the coefficients a, |3 ) and ( a, (3 ) can be written now in terms of the the real root 
yo[oc,^) of the polynomial 2v*+'^ + y"- — ]. For example, by following [[TT|, it is possible 
to write 



and 



: 1 - y^f (a^ « + 1 ) + PMa'o" - 1 )' - 2a|3 « - 1 : 



A/* -3 

+ 




(v«-l)(cx+(3-|3v«) ' 

where Vq = vo(a., |3). The root vo(cx:, |3) can be obtained as the value at x = 1 of a suitable 
analytic extension of the function defined within its convergence disk by the Taylor series 



y 2^ /-(pn+1)a-iV^p,^i ^ ^ (-2r /((a + |3)n + 1 - a) g'^ 



These extensions are finite sums of hypergeometric functions. 

The values of P(a, (3) are always positive integers. It is easy to show (see below) that 
there is always a minimum integer number PminfoC) (3) such that the spacing between 
consecutive allowed values of P(a, (3) is gcd(a, (3). Taking this fact into account and 
using vo((x, |3), |j.(a, |3) and ff((x, (3), it is possible to define a set of parameters that can 
be used to study the goodness of the approximations obtained by using the generalized 
peak counters introduced here. In particular, we will consider 



f n-. M-(a, |3)v/gcd(a, (3) 
o-(a, |3) 

amm(a, |3) = Pmm((x, |3)M-(a, (3) , 

amaxla, |3J = gcd(a, (3) — ^^^^^ — M.(a, |3J 



^ For the sake of simplicity, we will work without the projection constraint here because the leading terms 
for the mean and the variance are not sensitive to it. 
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and also 

log^ofa, (3) 



7rystairs(a, P) 



For our purposes, it suffices to consider those values of a and (3 for which the argument 
of the square root in 7rygauss('^? P) is positive. By using the following evident facts 

A/o(na,n|3) = yl'''^{oi, (3) , pmm(n(x, n|3) = npmin(a, |3) , 

|j,(n(x,n|3) = — |J.(a, |3) , cr^(na,n|3) = — cr^(a, |3) , 
n n 

it is possible to see that these parameters satisfy 

X(na,n|3) =x(a, |3) , amin(na,n(3) = Qminta, (3) , amax(na,n|3) = a^axta, (3) 

ystairs(Tl(X,n|3) = ystairs(a, |3) , ygauss(n(X, n(3) = ygauss((X, (3) , 

and, hence, they can be used consistently to assess the appropriateness of the different 
partitions. Although we will not give here a complete analytic proof, we provide enough 
numerical evidence to support the choice of P(3, 2) as the best peak counter (as was to 
be expected) in Figs. [8]and|9| where we show the values of x(oC) |3) and ystairsfo^) |3) for 
1 < a, (3 < 80. 



Several comments are in order now. First we want to comment on the role and meaning 
of amm('^> (3)- It is important to notice that all counters can be used to exactly reproduce 
the behavior of the entropy. This is so because they provide partitions of the set of black 
hole configurations. However, the behavior of the associated Gaussian approximations 
-available for all of them- differ for different choices for a, |3 and are not optimal in many 
instances. For example, before reaching the area value given by amm(a, |3), it is not true 
that the Gaussian approximation to the entropy can be understood as the sum of equally 
spaced Gaussian steps. In fact, the distance between consecutive steps in this regime is 
dictated by the values of p (non necessarily consecutive) that give non-zero values for 
[^''']G*(A^, 0; a, (3). When p > Pmin(c>c, |3) these values are separated by the gcd((x, (3) and 
are "as consecutive as possible". When gcd(a, |3) = 1 the value of Pminfc^) |3) corresponds 
to one plus the Frobenius number of the arithmetic sequence ak + (3 with k G N. This is 
given by [p^ (see also [20J) 

Pmin(a, (3) = 1 + a(a+ (3 - 1) . 

The value of amintcc, |3) by itself does not tell us anything about the quality of P( a, (3) as a 
counter. In addition to the threshold area, there is another relevant value Qmaxta, |3) -that 
can be roughly defined as the maximum area for which two consecutive steps can be dis- 
cerned in the gaussian approximation- that plays an important role. In fact, the interval 
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FIG. 8. Values of x(a, (3) for 1 < a, |3 < 80. The abscissa corresponds to |3. Darker colors represent 
the largest (better) values of xta, (3). The best choice is a = 3, (3 = 2 for which x(3, 2) ?a 35. 



length amax(«., (3) — amin(«,, (3) tells us how well the chosen Gaussian approximation cap- 
tures both the structure at the smallest area scales and the steps in the entropy. Though it 
is difficult to give a unique definition for amax(oc, |3), it is obvious that the broadening of 
the Gaussian steps signaled by the growth of the variance a^[a, (3) leads to the difficulty 
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FIG. 9. Values of Ttystairst^) P] for 1 < a, (3 < 80. Darker colors represent the largest (better) values 
of 7ty stairs ( CC) P ) • The best value corresponds to a = 3, |3 = 2 but other choices give values very 
close to 7tystairs(3,2). A similar picture can be drawn for 7tygauss(oi> P) that essentially coincides 
with this one for a > (3. When a < (3 the behavior changes because Ttygaussfa? P) becomes larger 
than 7ty or even complex. This does not change our conclusions about the goodness of P(3, 2). 



of separating two consecutive ones beyond 

amaxla, |3j = gcd[a, (3) — ^^i^^ — • 

This condition is derived by requiring that the width of a step is essentially equal to 
the distance from the previous one. Another quantitative criterion is provided by the 
exponents in ( V.8[ > and ( V.9[ > that tell us the decay rate of the staircase structure. The 
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inverse of the numerical coefficient of the area in the exponents of these expressions gives 
an order of magnitude estimate of amax(<^) |3). In practice the best choice of peak counter 
is the one giving both a low value of amin(c>c, (3) and a large value of a^axio^-, P)- The 
numerical evidence available tells us that the choice a = 3 and (3 = 2 is also optimal in 
this respect. 



VI.2. Peak substructures: Using two counters 



As we have discussed in the previous section, it is actually possible to partition the 
configuration space by using different types of linear counters P(a, (3). The best choice, 
as far as the description of the entropy structure is concerned, is P := 3K + 2N. The 
purpose of this section is to explore the effects of performing a further partition of the 
configurations corresponding to a single peak P = p by introducing an extra counter. The 
rationale behind this analysis can be immediately perceived by looking at the structure of 
a peak (see Fig. [lO]). As it can be seen, there are coherent substructures within each peak 
that are responsible for its asymmetric shape in a logarithmic plot. The most obvious and 
straightforward way to address the description of this substructure is to use an additional 
peak counter to further partition the space of configurations for a single peak. This can 
naturally be done by including extra variables in the generating functions. For example, 
when the projection constraint is ignored, the moment-generating function 

GJa/i,V2,s): 



1 -2Lr=i exp(-Sv/k(k + 2)) 



takes into account the contributions of the counters Pi = P(ai, (3i) and P2 = P(cx2, pi) 
This moment-generating function satisfies 



Ejexp(-sA) |pi,p2] 



and can be derived from the obvious generalization of the master generating function 
( |IL3| 



00 00 



G.(^'i,^2;xi,X2,...) := 1 -2}^}^vr''^^^Wr'^;^+P^x< . (VI.2) 



i=l n=l 



The black hole degeneracy spectrum for the configurations satisfying Pi = pi, P2 = pi 
and area a = qi i/pT + c|2a/P2 H is given by 



D*(a|pi,P2) = hV^T ■ ■ ■][A'i'^][vp]G*(^i,^'2;xi,X2, . . . ) . 

In Fig. |lO|we use Pi := P(3,2) = 3K + 2N and P2 := P(5,4) = 5K + 4N. This choice is 
favored by the arguments presented in |j8l and can also be understood by looking at the 
role of the Pell equation in the characterization of the area spectrum. A statistical descrip- 
tion of the subpeaks can be made by following the steps discussed in previous sections. 
In particular a Gaussian approximation can be obtained that, presumably, improves the 
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FIG. 10. Detail of the peak Pi = 3K+2N = 1 30. This logarithmic plot clearly shows the existence of 
several "subpeaks" labeled by the counter P2 = 5K+4N . Notice the asymmetric shape (mentioned 
in Fig. |2]). 
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one given by a single peak counter P(3, 2). However, as we do not expect the problems 
encountered above to be fixed by following this approach (in particular the underesti- 
mate of the value of y), we will not pursue it further. In any case, it is obvious that, by 
introducing extra counters and studying the statistical properties of the resulting peaks, 
one can get improved approximations for the entropy 

VII. CONCLUSIONS AND COMMENTS 

The main idea of the paper is to use the master generating functions, that encode 
the black hole degeneracy spectrum in an exact way, to derive statistical properties that 
can be used to describe and understand the detailed features observed in the black hole 
entropy. We have succeeded in reproducing, from a purely analytical point of view, the 
staircase of the entropy and its behavior as a function of the horizon area. In particular, 
we have shown how and why the steps disappear. The key element of our approach has 
been to use statistical properties of some subsets of black hole configurations (the "peaks" 
defined by P = 3K + 2N) to construct analytic approximations for the black hole entropy. 
It is very important to highlight the fact that we are not merely fitting the data but, rather, 
computing, a priori, the relevant statistical parameters of the peaks by employing the BH 
generating functions. In particular, we have given a procedure that can be used to get all 
the moments of the area distribution. 

It is also necessary to emphasize that the partitions of the space of BH configurations 
that we have exploited are exact, so there is no approximation introduced by the choice of 
the different "peak counters". The smoothing of the entropy consisting on adding Gaus- 
sian steps is a natural one that works remarkably well for small areas. Furthermore, there 
are important theorems in Combinatorics that guarantee the convergence of the individ- 
ual steps in the entropy (selected by the peak counter) to Gaussian distributions (after 
normalization). However, the sum of the Gaussian steps does not converge to the entropy 
because, as we have shown, the linear growth predicted the Gaussian approximation is 
slightly smaller than the actual one. It must be pointed out here that the numerical esti- 
mates for y, even in the crudest approximations, are remarkably good. In any case, the 
area range where the Gaussian approximation is reliable (that can be essentially obtained 
by comparing the actual growth given by the true y and ygauss) is large enough to trust the 
Gaussian approximation regarding the disappearance of the staircase structure (see Fig. 
|7|. Finally, we have discussed in Appendix |C] an alternate way of assessing the validity 
of our approach by looking at the pole structure of the Laplace transforms of the entropy 
and its Gaussian approximation. The comparison of both analytic structures tells us that 
they differ in their behavior far from the real axis. The numerical evidence available gives 
tantalizing evidence (encoded in an approximate periodicity that can be readily seen in 
Fig. [12| >, that prevents us from excluding a revival of the observed staircase structure for 
large area values. In any case, we deem this possibility quite unlikely. 

We have given numerical evidence to support the election of P = 3K + 2N as the best 
peak counter within the class P((x, (3) = aK + (3N, a, (3 G N. Its usefulness is justified 
by the fact that in the low area regime the variance of the distributions associated with 
the peaks is much smaller than the separation of the mean areas corresponding to con- 
secutive peaks. This explains why a staircase structure must be seen in this regime. The 
actual comparison of the exact entropy values and the prediction given by our model is 
very compelling. In any case, it is conceivable that other functions, more general than 
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the counters P ( a, (3 ) that we have discussed, can be defined in order to better understand 
and approximate the large area behavior of the entropy and the value of y. 

Although most of the computations carried out in the paper have made use of the DL 
prescription to obtain the black hole entropy [12J, our methods can be extended in a com- 
pletely straightforward way to other countings such as the SU(2) proposal of BSllll. The 
numerical details differ from the ones that we have presented above but the qualitative 
conclusions -which are independent of the concrete form of the projection constraint or 
equivalent conditions- remain unchanged. In particular, when the condition that plays 
the role of the projection constraint in this case is ignored, the moment-generating func- 
tion is lHOl 



1 - ^ + 1 )^"'+' exp(-s v/k(k + 2)) 

By following the methods described in the paper step by step, one finds that yo = 
[0.73996900 ■ ■ ■) is now the smallest real root of the polynomial + — 2y^ — 2y^ + 1 . 
The values for \x and cr^ are obtained from the function 



^,s) = 1 -2^(k+l)v3'^+^exp(-sv/lc(k + 2)) 



xENP 

k=l 



They are p. = (0.34980945 ■ ■ ■ ) and = (0.00010926 • • ■ ). The staircase and Gaussian 
approximations to the entropy lead to the following values for the Immrizi parameter: 
Try stairs = [0.86088860 ■ ■ ■) and Tty gauss = (0.86100438 ■ ■ ■ )• These have to be compared 
with TiyENP = (0.86100642 ■■• ) obtained in |21||. The analytic structure of the Gaus- 



sian approximation is similar to the one shown in Fig. 11 and, then, the behavior of 
the entropy is essentially the same as in the DL case. It is obtained from this one by 
changing the values of the parameters yo, and cr^. The projection constraint can be in- 
troduced in the same way as before and, as expected, only configurations with even P 
contribute in this case. The relevant singularities in the Laplace transform are located at 
s,± = (0.81 057903 • • • ) T ( 1 7.96628553 ■ ■ ■ )i and s± = [0.84839767 ■ ■ ■ ) T (8.98324891 ■ ■ ■ )i. 
This leads to a doubling of the size of the steps and the persistence of the staircase struc- 
ture for larger values of the horizon area in the case when the projection constraint is 
taken into account. 

A last comment that we want to add is related to the values of the areas for which the 
steps in the entropy and the peaks in the BH degeneracy spectrum cease to be seen. As 
the entropy is obtained by integrating the degeneracy spectrum and taking the logarithm 
of the resulting sum, it is to be expected that the effective disappearance of the steps takes 
place for smaller values of the areas than the disappearance of the corresponding peaks in 
the degeneracy spectrum. This must be taken into account in order to correctly interpret 
the meaning of the substructures found in the behavior of the entropy as a function of the 
area. 
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Appendix A: Computation of the moments when the projection constraint is not considered 



In this Appendix we will derive the mean and variance of the area distribution in the 
non-projection constraint setting by using (IIL8). They are given by 




2y5 + ^3 



and (T^p = [A^ | p] - ^Jp, where 

/ 



EjA^lp] 



1 



V 



(2^/5 + A/3 -1)3 



Here, as in the main body of the paper, we use the notation 



a* 



2^5 + ^3 _ 1 



The computation of the coefficients in the Taylor expansions about y — that appear in 
the previous formulas can be easily carried out, for instance, by using the Cauchy integral 
theorem 



[ym^] 



1 

27ri 



() 



d^. 



where C is an index-one curve surrounding the origin y = (and leaving the remain- 
ing singularities of the integrand outside). The pole at the origin has order p + 1 so, in 
practice, it is better to compute the integral by moving the contour radially outwards and 
picking up the contributions of the remaining singularities of f ( a/) / y^^^ . This is useful be- 
cause they are, in many cases, poles of a fixed, p-independent order with a p-dependent 
residue. The value of a*p can be easily obtained by using this procedure: 



1 

2m 



() 



1 



1 



= - ^ Res 



^ A/p+i 2y^ + v3 - 1 

1 



dy 



(2^5 + ^3 _ 1 



L 

i=-2 



ly} 



^Qy\-r^)y\ 



where the yi are the five different roots of the polynomial ly'^ + "v^ — 1 . We have used 
the convention that y^ and v i are complex conjugate of each other. The single real root, 
y^i — [0.77039825 ■ ■ ■ )/ is the one with the smallest module (a closed expression for y^ in 
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terms of hypergeometric functions is given in IITOll '). The previous expression shows that 
the asymptotic behavior for large values of p of (X*p is given by 



2yi 1 

a*p ~ ^ T — , p ^ oo. 



There are other alternative expressions for the coefficients a*p. For example they satisfy 
the simple recursion relation 

o(-*'P+5 = oc*p+2 + 2a*p, p > 0, a*o = 1 , cx:*i = a*2 = oc^s — oi^4 — . 
They can also be written in terms of binomial coefficients as 

The numerators appearing in the expressions for the mean and the variance can also be 
obtained in closed form by using the Cauchy integral formula. However the expressions 
thus obtained are not very illuminating so, here we will only give their leading asymp- 
totic behavior for large values of p. In both cases this is given by a term linear in p: 

^,p ~ ^ ■ p + Ho = (0.34959022 ■■■)■?+ {0.000^ 9724 ■ ■ ■ ) 
o^p-ct2-P = (0.00009817---)-P, 

where \l, [Xq, and are given in terms of yo: 



.k=l 




where 
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f Syl{20yl + 3] 2(5v| + l) \ ^ ^ 
= (3 - 2|a. • (Xo , 



/ 12(20vg + 3)(1-vg)^ 6[}-yin5yl + }, . /UnTTITv 




k=l 



'0 ». "0 

= (0.00023608 




The results obtained by this method are completely equivalent to those presented in Sec- 
tion |Vj Higher moments can be computed by using the same procedure described here. 
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Appendix B: Computation of the moments when the projection constraint is considered 



As we have discussed in subsection IIL2 the incorporation of the projection con 



straint requires us to introduce of an extra variable z in the moment-generating function 
G{y,s]z]. To compute the moments in the area distribution we have find out some s- 
derivative, evaluated at s = 0, of G[y, s; z). This derivative gives rise to a new function 
f[y.)Z) whose coefficient fp = [z°][A''']f(z, v) allows us to determine the moments in the 
area distribution. This process defines a function f[y) = Y.p fp^^ of the single variable y 
that is desirable to have in closed form. We describe in this subsection how this can be 
done. What we need to find is 



p=0 



dz. 



whenever f ("v, z) is an analytic function of z for all the values of "v in a neighborhood 
of "v = 0. If [v'llf (z, "v) is a Laurent polynomial at z = we can choose any contour 
surrounding z = for Ci (as long as it is piecewise smooth and of index one). Now if 
f ("v, z) is analytic at = for all the values of z G Ci the previous expression can be 
rewritten as 



p=0 



2m 



() 



dz 1 



2m 



() 



C2 



dy 



where the contour Cz (again piecewise smooth and of index one) is now chosen in such 
a way that, for each z G Ci the only singularity surrounded by it is v = 0. We can now 
exchange the order in the integrals to get 



y{y) 



p=0 



2m 



() 



dy 1 

A/P+i 2m 



() 



dz 



f{y,z] 



1 

2m 



() 



dz 



f^,z) 



(B.l) 



whenever the last integral is an analytic function of a' in an open neighborhood of = 0. 
If the integration contours are chosen according to the prescription described here, there 
are (generically) singularities of z i-^ f[y,z)/z inside Ci (in most of the cases z = and, 
eventually, others coming from f(v, z); which will be functions of y). The residues at 
these singularities give us a closed-form expression for 



Ffvl = 



1 

2m 



() 



dz 



f(v,z) 



In many practical situations it is convenient to choose a unit, positively oriented, circum- 
ference for Ci. This is so because, for this choice, it is possible to simplify expressions of 
the type 
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2m 



() 



dz 



z'^ + z '^)g(v,z + z 



1 

m 



() 



dz 



z''g(A/,z + z" 



kG N, 



as can be easily seen by performing the change of variable z i-^ z^ in the second term 
of the first integral. In the particular case of interest, by setting s = in the master 
generating function G [y, s; z), we find 



f(A/,z) = G^,0;z) 



1 



— zy^ — z + y^ 



1-L 



oo 
k=l 



z^ + z-^]y^^+^ 



zJ-y^ + z^y^ — Izy^ — zy^ — z + y^ + y' 



3 • 
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In order to apply the previous procedure we first take a unit, positively oriented, circum- 
ference for Ci . With this choice, and once a suitable C2 contour is picked, the two relevant 
poles in the integrand (B.l I are z — and 

_^+2y^ + y^- y/{y^ - 1 ) (a/4 + v2 + 1 ) (2^5 + [2y^ + y^ + ]) 



2y^ (y^ + r 



so that 



1 +y^ 



1 + V 



1 



2a/5 + A/3 - 1 / V 2a/5 + A/3 + 1 



A/3+1 



Some facts are evident at this point. In particular, the coefficient ccp := [A/i'lFfA/) is zero for 
odd values of p. For large even values of p the asymptotic behavior of cXp is controlled by 
the singularities y = ±a/o and is 



1 



1 



CX2q 



The mean yi^ is computed by following the same steps as in the non-projection con- 
straint case. In particular (Xp ■ [ip is given by the [y^] coefficient of 



2A'^fA'^ 



1 



;i+A/^FRfA/; 



2a/^(1 + a/2) 



k=l 



QiM + v/RM 

2y^n +A/2) 



y] + JR[y] 



where 



R[y) 



(a'2 - 1 ) (a/4 + A/2 + 1 ) (2^5 ^ A/3 + 1 ) (2a/5 + A/3 

2y^ + y^ + ^ 
2 + A/2 - 2a/^ - 7y^ - 6y 
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For large even values of p, the asymptotic behavior of [i^, is also controlled by the singu- 
larities y = ±A/o. In particular, by using the identities 



2v^(l+A/2 



^2A/8 + A/g + 1 + v/R(±^ 
it is straightforward to show that 

fX2q ~ 2|j,q , q 
where, as in the non-projection constraint case. 



4A/2 



2a/« + v^ + 1 



00 . 



IOa/2 + 3 



2^ Vk(k + 2)^ 



3k+2 




,3k 
'0 > 



,k=l 
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Finally, (cTp + |j.p) ■ ocp is given by the [y'P] coefficient of 



CO k] —1 



+ 

y) VR(v) [Qs[y) ./RM - Qi [ymy] 

ki-k2 



+ 



1 + A/2 



(^Qjiy) + 3(k2 - ki ) VQ2(^) x/RM - (k2 - ki ) Vr(a/; 



ki +k2 



1 



(QzM -3(ki +k2)VQ2Mv^- (ki +k2)VR(V 



where 

Q3(^) 
Q4(V) 

Q5(V) 



1_^2^^4_3^6_3^8_5^10^3^12^3^14^-,1^16^7^18^7^20_^22_4^26 

3+3y^W-9y^-27y^-27y^°+5y^^+47y^^+8]y^^+65y^^+5y^°-55y^^-64y^^-28y^'^ 
-2(3+3a'^+a/4-6a/^-24a'^-26a/^°-a/^^+23a/^V55a'^^+64y^^+28a'^°) . 



The asymptotic behavior of 02q can be obtained in a straightforward (albeit tedious) way 
from the preceding expressions. 



Appendix C: Analytical properties of the Laplace transform of the entropy in the Gaussian 

approximation 

In order to understand the main features of the Gaussian approximation for the en- 
tropy, in particular its asymptotic behavior as a function of the area, it is very convenient 
to rely on Laplace transform methods. This is so because the asymptotics of a function 
can be understood in many cases by looking at the analytic structure its Laplace trans- 
form, specifically the locations of its singularities and their type. In our case, and despite 
the fact that the behavior of the Gaussian approximation can be roughly understood from 
the features of the individual steps, the detailed behavior of the sum is much harder to 
get. In particular the realization of the fact that the value of the Immirzi parameter is not 
exactly recovered in this approximation really demands a detailed analysis for which the 
use of Laplace transform methods is especially appropriate. The arguments given below 
refer to P(3, 2) but can be trivially extended for other peak counters P(a, (3), defined in 
Section |Vl| as long as 

1 + ?^logv„(«,p)>0. 

^2((X, (3) 
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s — plane 



FIG. 11. We represent here the regions relevant to the discussion of the convergence of the func- 
tional series ( C.2| . We also plot the position of the relevant singularities of the Laplace transform 
(there is an infinite number of other isolated singularities also located on the hyperbola which are 
not relevant for the analysis presented in the paper). 



Let us consider first the non-projection constraint case 



£(expSr",s) 



1+ -0 flfl-erff^ 
s + (i0v^ + 3)s^v^ ^ [via 
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1 +erf ( — T=-^^ 



Via 



The well known asymptotic behavior of the error function 



exp 



(C.l) 



M-s p 



erf (x) ^ 1 



7tX 



X > 1 , 



guarantees the convergence of the first series in ( C.l[ >, hence, the first to terms have a very 
simple singularity structure: just a simple pole at s = 0. Let us concentrate then in the 
last series in dCTl) 



^ 1 + erf — ^-^Vv 



exp 



0^ 



(C.2) 
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FIG. 12. Comparison between the absolutes values of the Laplace transforms £(exp S*, Tty + iy) 
(continuous line) and L[exp ST^^^, Tty + iy) (dashed line). 



By using now the following asymptotic formula for the error function for complex values 
of its argument 



1 / °° 

•■[^-er^[z]]-^{^ + y[-^y 



(2m- 1)!! 

(2z- 



2-\m 



, z^oo, |arg(z)|< 



371 

T 



it is possible to prove the point convergence of the series (C.2) in the wedge | arg(s — 
\x/ a^) I < 37t/4. In the region | arg(s — yi/c^] \ > 3n/4 the asymptotic behavior is given by 



erf(z) 



-1 



37t 



, z^ oo |arg(z)| > 

TTZ 4 



In order to have point convergence in this case for (C.2) we have to demand ( V.4[ ) 



Re 



2 

(yS^- M-S-logVo) =0. 



The divergence of the Laplace transform on the hyperbola defined by this condition can 
be absorbed in the singularities of the function 



^ Vexp (p( 

p=0 ^ 
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1 — exp ^^y^ — |J.S — logA^o 
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in the region | arg(s — |J./o"^)| > 37r/4 (see Fig. 11 ). These singularities control the asymp- 
totics of the Gaussian approximation to the entropy as discussed in Section |V| When the 
projection constraint is taken into account it is straightforward to identify the function 
Fging that encodes the singularities of the Laplace transform of the Gaussian approxima- 
tion to the entropy 



T'smg(s) = Y_ exp \^2q (^-^ ^ (^^P {^^^ ~ -^^^ - 21og a/q 



We end this appendix by giving an additional way to compare the exact behavior of the 
entropy and the Gaussian approximation. The Laplace transform inversion formula tells 
us that given the Laplace transform of a function £ (f ; s) it is possible to recover f by using 
the inversion formula 

^(f-.sje'^Ms. 

J xq— i-oo 



f(a) = 

Inx 



The integration contour can be taken to be a line Re(s) = xq, parallel to the imaginary 
axis leaving all the singularities of the integrand to its left. The Laplace transform of the 
exact entropy is 

/L(expS*,s' 



l-2Lr=i exp(-sv/k(k + 2)) 



If we now compare the Laplace transforms £(exp S*, Try + iij) and £(exp S^"^^, Try + iij) 



by plotting their absolute values as functions of y we get the result shown in Fig. 12 
As it can be seen there is a remarkable agreement between both plots for small values of 
y, however this agreement disappears for larger values. This is due to the fact that the 
real parts of the poles of the Laplace transform of the Gaussian approximation do not 
accumulate because they are located on hyperbolas as explained in this Appendix and in 
Section rvl 
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